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Abstract 

We  study  the  heavy  traffic  regime  of  a  multiplexer  driven  by  correlated  in¬ 
puts,  namely  the  M\GI\oo  input  processes  of  Cox.  We  distinguish  between 
M\GI\oo  processes  exhibiting  short  or  long-range  dependence,  identifying  for 
each  case  the  appropriate  heavy  traffic  scaling  that  results  in  non-degenerate 
limits.  As  expected,  the  limits  we  obtain  for  short-range  dependent  inputs 
involve  the  standard  Brownian  motion.  Of  particular  interest  though  are  our 
conclusions  for  the  long-range  dependent  case:  The  normalized  queue  length 
can  be  expressed  as  a  function  not  of  a  fractional  Brownian  motion,  but  of 
some  other  stable  non-Gaussian  self-similar  process.  Thus,  the  M\GI\oo  pro¬ 
cesses  serve  as  an  example  demonstrating  that,  within  long-range  dependence, 
fractional  Brownian  motion  does  not  assume  the  ubiquitous  role  that  its  coun¬ 
terpart,  standard  Brownian  motion,  plays  in  the  short-range  dependence  setup, 
and  that  modeling  possibilities  attracted  to  non-Gaussian  limits  are  not  so  hard 
to  come  by. 
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1  Introduction 


The  apparent  presence  of  long-range  dependence  and  self-similarity  in  network  traf¬ 
fic,  that  has  been  established  in  numerous  studies  of  WAN  [21],  Ethernet  [13]  and 
VBR  video  [3]  measurements,  raises  the  need  to  revisit  various  performance  anal¬ 
ysis  and  design  issues  in  this  new  regime.  As  some  recent  experimental  work  ([7]) 
already  suggests,  long-range  dependence  has  a  tangible  and  adverse  effect  on  queue¬ 
ing  measures  such  as  buffer  overflow  probabilities,  and  consequently  its  presence  is 
not  to  be  overlooked  or  underestimated. 

Roughly  speaking,  this  phenomenon  of  long-range  dependence  amounts  to  cor¬ 
relations  in  the  packet  stream  spanning  over  multiple  time  scales,  which,  although 
individually  rather  small,  decay  so  slowly  that  they  are  non-summable  in  the  long 
run,  thus  affecting  performance  in  a  drastically  different  manner  than  that  predicted 
by  traditional,  summable  correlations,  Markovian  models.  In  an  effort  to  gain  some 
understanding  into  long-range  dependence,  and  in  particular  into  the  fundamentals 
of  its  impact  on  queueing,  various  traffic  models  that  account  for  non-summable 
correlation  patterns  have  been  proposed.  These  models  include,  among  others,  frac¬ 
tional  Gaussian  noise  inputs  [1],  fractional  Brownian  motion  [16]  and  on/off  sources 
with  Pareto  interrenewal  periods  [5],  with  early  results  yielding  buffer  asymptotics 
of  radically  different  nature  than  the  exponential  tails,  that  typically  characterize 
Markovian  models. 

In  this  paper,  we  consider  the  class  of  M\GI\oo  input  processes,  which  fall  in 
the  category  of  models  that  can  account  both  for  short  and  long-range  dependent 
behaviour.  An  M\GI\oo  input  process  is  understood  as  the  busy  server  process 
of  a  discrete-time  infinite  server  system  fed  by  a  discrete-time  Poisson  process  of 
rate  A  (customers/slot)  and  with  generic  service  time  a.  Such  M\GI\oo  processes 
have  already  been  used  by  Paxson  and  Floyd  to  successfully  model  WAN  traffic  in 
[21].  What  is  more,  a  very  attractive  justification  for  M\GI\oo  modeling  is  given 
by  the  fact  that,  as  it  is  argued  in  [14]  and  [20],  these  processes  emerge  simply 
by  combining  traffic  generated  by  on/off  sources  with  a  Pareto  distributed  activity 
period,  and  letting  the  number  of  sources  go  to  infinity,  so  as  to  obtain  a  behaviour 
identical  to  that  of  a  M\GI\oo  input  stream  with  Pareto  distributed  a. 

The  M\GI\oo  processes  are  very  flexible  and  allow  for  direct  control  of  their  cor¬ 
relation  patterns  through  a  (see  Proposition  2.1  and  (2.1)).  They  are  also  extremely 
tractable;  a  variety  of  buffer  tail  probabilities  for  a  multiplexer  fed  by  M\GI\oo  in¬ 
puts  are  obtained  in  [18]  and  [20].  Surprisingly  enough  though,  in  the  long-range 
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dependence  setup,  these  results  provide  hyperbolic  buffer  asymptotics,  altogether 
different  from  the  Weibullian  tails  obtained  under  the  fractional  Brownian  motion 
assumption.  Here,  we  will  try  to  further  explore  this  fact,  this  time  by  devel¬ 
oping  some  queueing  theory  for  an  M\GT\oo  process  feeding  a  multiplexer  that 
operates  in  heavy  traffic.  Our  motivation  for  going  to  heavy  traffic  is  twofold  : 
First,  given  that  under  short-range  dependence  heavy  traffic  analysis  has  offered 
useful  standard  Brownian  motion  characterizations  of  queueing  networks,  it  would 
be  interesting  to  try  to  extend  the  theory  to  the  long-range  dependence  setup.  This 
would  possibly  furnish  evidence  demonstrating  that  fractional  Brownian  motion,  the 
long-range  dependent  analog  of  standard  Brownian  motion,  plays  a  similar  central 
role  in  the  modeling  of  long-range  dependent  traffic.  Secondly,  for  the  particular 
case  of  M\GI\oo  inputs,  a  heavy  traffic  analysis  would  help  elucidate  the  observed 
differences  in  buffer  asymptotics  between  M\GI\oo  and  fractional  Gaussian  noise 
inputs.  Unless  these  differences  are  due  to  some  fundamental  structural  discrepan¬ 
cies,  one  would  expect  that,  in  the  heavy  traffic  regime,  both  models  would  collapse 
to  a  fractional  Brownian  motion  characterization,  similarly  to  what  happens  under 
short-range  dependence,  where  different  models  eventually  collapse  to  a  description 
involving  the  standard  Brownian  motion.  On  the  other  hand,  if  the  abovementioned 
differences  are  due  to  some  different  salient  features  of  the  M\GI\oo  and  fractional 
Gaussian  noise  process,  it  would  be  reasonable  to  expect  that,  despite  the  asymptot¬ 
ically  identical  correlation  patterns,  the  differences  would  carry  over  and  manifest 
themselves  even  more  clearly  in  the  heavy  traffic  regime. 

The  results  of  our  investigation  confirm  the  latter  statement.  First,  we  find 
that,  under  short-range  dependence,  the  class  of  M\GI\oo  inputs  belongs  to  the 
domain  of  attraction  of  the  standard  Brownian  motion,  as  expected.  Next,  and 
most  importantly,  we  show  that  under  long-range  dependence,  in  particular  for  a  a 
in  the  domain  of  attraction  of  a  non-normal  stable  law,  the  M\GI\oo  process  is  not 
attracted  to  a  fractional  Brownian  motion.  Instead,  the  normalized  heavy  traffic 
queue  length  is  expressed  as  a  function  of  a  non-Gaussian  a-stable  self-similar 
process.  This  result  underlines  the  fundamentally  different  nature  of  the  long-range 
dependent  M\GI\oo  processes  and  also  points  to  the  fact  that  fractional  Brownian 
motion  does  not  necessarily  play  the  same  key  role  that  standard  Brownian  motion 
assumes  under  short-range  dependence.  Within  long-range  dependence,  there  seems 
to  be  a  choice  for  distinct  modeling  possibilities;  and,  as  our  study  suggests,  it  is 
not  at  all  difficult  to  encounter  a  rather  simple,  potentially  useful  model,  that  is 
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attracted  to  non-Gaussian  limits. 

The  rest  of  the  paper  is  organized  as  follows:  The  class  of  M\GI\oo  input  pro¬ 
cesses,  along  with  some  of  their  properties,  is  introduced  in  Section  2.  We  present 
our  multiplexer  model  in  the  heavy  traffic  regime  in  Section  3.  Our  heavy  traffic 
results  are  then  stated  in  Section  4,  while  Section  5  addresses  their  implications. 
Finally,  we  give  an  outline  of  our  proofs  in  Section  6. 

A  few  words  about  the  notation  used  here.  All  RVs  are  defined  on  some  prob¬ 
ability  triple  (O,  T,  P),  with  E  denoting  the  corresponding  expectation  operator. 
Weak  convergence  is  denoted  by  =>-.  We  write  /  ~  g  when  lim  f(x)/g(x)  =  1,  and 

fr  x  gT  when  lim  fT(x )  —  lim  g(x),  x  G  IR. 
r—¥  oo  r— >-oo 

2  M\GI\oq  processes 

As  we  intend  to  model  short  and  long-range  dependence  through  the  busy  server 
process  {bn,  n  =  0, 1, . . .}  of  a  discrete-time  M\GI\oo  system,  we  collect  in  this 
section  some  pertinent  facts  about  { bn ,  n  =  0, 1, . . .},  that  have  already  been  estab¬ 
lished  elsewhere  ([6],  [18],  [19]),  and  will  be  used  in  the  sequel. 

To  describe  the  process  {6n,  n  =  0, 1, . . .},  consider  a  system  with  infinitely  many 
servers  :  During  time  slot  [n,n  +  1),  fJn+i  new  customers  arrive  into  the  system. 
Customer  j,  j  —  1, . . . ,  /3n+ 1,  is  presented  to  its  own  server  and  begins  service  by  the 
start  of  slot  [n  + 1,  n  +  2);  its  service  time  has  duration  crn+i,j-  Then  bn  denotes  the 
number  of  busy  servers,  or  equivalently  of  customers  still  present  in  the  system,  at 
the  beginning  of  slot  [n,  n  +  1),  with  b  denoting  the  number  of  busy  servers  initially 
present  in  the  system  at  n  =  0. 

The  .CV-valued  RVs  6,  {/?n+i,  n  =  0, 1, . . .}  and  {crn,j,  n  =  0, 1, . . . ;  j  =  0, 1, . .  .} 
satisfy  the  following  assumptions:  (i)  The  RVs  are  mutually  independent;  (ii)  The 
RVs  {/3n+i>  n  =  0, 1,...}  are  i.i.d.  Poisson  RVs  with  parameter  A  >  0;  (iii)  The 
RVs  {cr nj,  n  —  1, . . . ;  j  =  1, 2, . . .}  are  i.i.d.  with  common  pmf  G  on  {1, 2, . . .}.  We 
denote  by  cr  a  generic  W-valued  RV  distributed  according  to  the  pmf  G.  Throughout 
we  shall  assume  that  E  [cr]  <  oo. 

It  is  also  helpful  to  introduce  the  forward  recurrence  RV  a  associated  with  cr. 
This  follows  the  equilibrium  distribution  of  cr,  i.e., 

^  00 

Tn  =  P  [5  >  n]  =  eToT  ^  P  [cr  >  j] ,  n  =  0, 1, . . .  (2.1) 

L  J  j=n 
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and,  as  part  2  of  the  following  proposition  shows,  it  directly  controls  the  correlations 
in  the  {bn,  n  =  0, 1, . . .}  stream. 

Proposition  2.1  The  busy  server  process  {bn,  n  —  0,1,...}  admits  a  stationary 
ergodic  version,  still  denoted  by  {bn,  n  =  0, 1, . . with  the  following  properties: 

1.  For  each  n  =  0, 1, . . .,  the  RV  bn  is  a  Poisson  RV  with  parameter  AE  [a], 

2.  Its  covariance  function  is 


co v(6«,  bj)  =  AE  [a]  r^j,  n,j  =  0, 1, . . . 


where  Vn,  n  =  0, 1, . . .,  is  given  by  (2.1)  . 

Additional  information  about  the  stationary  busy  server  process  {bn,  n  =  0, 1, . . .} 
is  contained  in  the  characteristic  function  of  the  partial  sums  of  the  process,  i.e.,  in 

f  /  r  \l 


$r(0)  =  E 


exp 


id  Y1  bn 


r  =  1,2, ... ,  9  £  1R 


L  \  n=i  / J 

Since  this  quantity  plays  instrumental  role  in  our  results,  we  will  list  below  the 
complete  expression,  whose  derivation  can  be  found  in  [18]  and  [19]  : 

Proposition  2.2  For  each  r  =  1,2,...,  and  9  £  1R,  it  holds  that 

ln$r(0)  =  -AE  [<r]  (l  -  E  j  +  AE  J(r  _  CT)+e^j 

— Ar  +  A  (l  -  e-^y1  E  -  l]  (2.2) 

Having  presented  the  necessary  facts  about  the  M\GI\oo  busy  server  process,  we 
next  describe  our  heavy  traffic  setup. 


3  The  multiplexer  in  heavy  traffic 

We  view  the  multiplexer  as  a  discrete-time  single  server  queue  with  infinite  buffer 
capacity,  operating  at  a  constant  rate  and  in  a  first-come  first-served  manner:  Let 
Qn  denote  the  number  of  cells  remaining  in  the  buffer  by  the  end  of  slot  [n  —  l,n), 
and  let  bn+i  denote  the  number  of  new  cells  which  arrive  at  the  start  of  time  slot 
[n,  n  +  1).  If  the  multiplexer  output  link  can  transmit  c  cells/slot,  then  the  buffer 
content  sequence  { Qn ,  n  =  0, 1, . . .}  evolves  according  to  the  Lindley  recursion 

Qo  =  0;  Qn+i  —  \Qn  +  frn+l  —  c]+,  n  —  0, 1, . . .  (3-1) 
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The  stationary  busy  server  process  of  a  discrete-time  M\GI\oo  system  we  dis¬ 
cussed  in  section  2,  will  now  be  used  to  model  the  arrival  process  { bn ,  n  —  1,2,...,...}, 
thus  providing  us  with  a  very  simple  way  of  accounting  for  time  dependencies  in 
the  cell  input  stream.  So,  with  the  M\GI\oo  process  { bn ,  n  =  1,2, be¬ 
ing  characterized  by  the  pair  (A,  a),  the  average  input  rate  to  the  multiplexer  is 
E  [bn]  =  AE  [cr] ,  and  it  can  be  shown  that  the  system  is  stable  if  AE  [a]  <  c. 

Our  goal  here  is  to  examine  the  behaviour  of  a  system  that  is  almost  fully  utilized. 
This  typically  involves  obtaining  limiting  expressions  of  properly  rescaled  quantities 
of  interest,  as  the  traffic  intensity  p  tends  towards  its  critical  value  1.  Although, 
in  general,  it  may  be  possible  to  derive  different  flavors  of  heavy  traffic  results,  one 
rather  standard  technique  for  studying  the  behaviour  of  a  queue  in  heavy  traffic 
is  to  consider,  instead  of  a  single  queueing  system,  a  family  of  queueing  systems, 
indexed  by  an  integer  parameter  r  =  1,2,  —  By  taking,  for  the  rth  system,  the 
input  and  the  service  processes  to  suitably  depend  on  the  parameter  r,  it  is  possible 
to  ensure  that,  as  r  —>•  oo,  the  corresponding  traffic  intensity  pr  goes  to  1.  Adapting 
this  approach  to  our  context,  we  consider,  instead  of  just  (3.1),  a  family  of  Lindley 
recrusions  where,  for  each  r  =  1,2,... 

QS  =  0;  Qrn+i  =  [Qn  +  K+i-c}+ ,  n  =  0,l,...  (3.2) 

Note  that  the  stationary  M\GI\oo  busy  server  processes  { brn ,  n  =  1, 2,  •  •  -},  r  = 
1,2,...,  above  are  described  by  the  Poisson  arrival  rate  and  service  RV  pair  (Ar,  a). 
That  is,  in  our  parametrization,  every  member  of  this  family  of  queueing  systems 
depends  on  r  only  through  the  input  process,  and  in  particular  its  Poisson  rate 
parameter  Ar,  while  the  service  RV  a  and  the  multiplexer  release  rate  c  remain  the 
same. 

It  can  be  shown  that  the  Lindley  recursions  (3.2)  admit  the  following  equivalent 
representation,  which  is  useful  for  establishing  heavy  traffic  limit  theorems:  By 
setting,  for  each  r  =  1, 2, . . ., 

Sr0=  0;  =  n  =  l,2,... 

i=i 

we  can  write 

Qn  =st  S„  —  inf  {5} ,  j  =  0,l,...,n},  n  =  l,2, ... 

Before  putting  ourselves  in  the  heavy  traffic  regime  by  letting  the  traffic  intensity 
go  to  1,  it  is  necessary  to  rescale  in  a  meaningful  manner,  both  in  time  and  in  state 
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space,  the  processes  of  interest.  We  do  so  by  introducing  the  scaling  £r,  r  =  1, 2, . . 
to  define  the  continuous  time  processes 


r  _  %r] 

(ft  =  -M: 

Sr 


t  >  0 


and 


so  that 


t  >  0 


ft = *1  - t-°  (3-3) 

for  every  r  =  1,2,  —  The  stationary  version  can  be  obtained  by  letting  t  -»  oo  and 
then  r  — >  oo,  in  order  to  drive  the  system  to  heavy  traffic.  Here,  we  will  interchange 
the  order  of  the  limits,  taking  r  — >  oo  first  and  t  — t  oo  afterwards,  to  again  collect 
the  same  stationary  RV  (see  [11]).  For  a  classical  GI\GI\1  queueing  setup,  where 
second  moments  are  always  assumed  to  be  finite,  the  scaling  that  gives  rises  to 
non-degenerate  limits  is  Cr  —  sjr.  In  our  M\GI\oo  situation  though,  that  allows  for 
an  infinite  second  moment  of  <r,  the  question  arises  as  to  what  should  the  suitable 
scaling  be  when  the  input  process  exhibits  long-range  dependence.  Although  the 
answer  is  not  immediately  obvious,  it  is  easy  to  see  that,  in  all  cases,  to  avoid 
collecting  only  a  law  of  large  numbers  result,  any  candidate  Cr  should  obey  the 
following  necessary  condition  : 


Condition  3.1  The  sequence  £r,  r  =  1, 2, . . satisfies 

Cr 

lim  Cr  —  +oo  and  lim  —  =  0. 

r— >oo  i — >-oo  t 

Without  any  further  requirements  on  £r,  we  may  now  enforce  the  heavy  traffic 
assumption  below,  which  guarantees  that,  as  r  -4  oo,  the  family  of  queueing  systems 
described  by  (3.2)  gradually  approaches  instability. 

Assumption  3.1  The  sequence  Xr,  r  =  1, 2, . . .,  satisfies 

T 

lim  (ArE  [cr]  —  c)  —  =  —7 

r— kx>  J  Cr 

for  some  7  >  0. 
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4  The  heavy  traffic  results 


In  the  heavy  traffic  setup  of  the  previous  section,  our  first  task  is  to  identify  the 
appropriate  scalings  £r,  r  =  1,2,...,  for  the  various  choices  of  the  distribution  of  the 
RV  cr,  that  controls  the  correlations  in  the  input  cell  stream.  Clearly,  our  aim  is  to 
obtain,  as  r  — >  oo,  a  non-trivial  limiting  process  (s^0,  t  >  0}.  Fully  characterizing 
this  process,  especially  in  the  somewhat  intriguing  case  of  long-range  dependence, 
will  then  provide  information  about  the  buffer  content  distribution  through  (3.3). 

With  Assumption  3.1  being  enforced  throughout,  we  distinguish  between  three 
different  cases,  which  correspond  to  different  tail  behaviour  of  the  distribution  of  a. 
It  is  worth  pointing  out  here  that,  for  each  case,  the  scaling  ( r  is  essentially  unique, 
in  the  sense  that  any  other  scaling  that  yields  a  non-degenerate  limit  is  related 
to  Cr  according  to  ('r  =  0((r). 

Our  first  result  addresses  the  case  where  the  M\GI\oo  process  is  short-range 
dependent,  i.e.,  Tn  <  +oo,  or,  equivalenty  (see  [18]),  E  [cr2]  <  +oo. 

Proposition  4.1  (Short-range  dependence)  If  E  [cr2]  <  +oo  then  with  (r  —  y/r, 
r  =  1, 2, . . it  holds  that 


{srt,t>  0}  {—7 1  + 

where  {Bt,t  >  0}  is  a  standard  Brownian  motion 


1 


C(1  +  2  £  rn)  Bu  t  >  0} 

n= 1 


The  next  result  deals  with  what  turns  out  to  be  a  cross-over  case  between  short 
and  long-range  dependence.  Despite  the  fact  that  E  [a2]  =  +oo  and,  as  a  con¬ 
sequence,  the  sum  of  the  correlations  is  Fn  =  +oo,  a  is  still  attracted  to  a 

Gaussian  RV,  while  X)n=i  is  a  slowly  varying  function. 

Proposition  4.2  Let  E  [cr2]  =  +oo  and  E  [1  [cr  <  x]  cr2]  ~  2E  [cr]  h(x),  where  h(x) 
is  a  slowly  varying  function.  Then  with  (r  =  y/rh{r),  r  =  1,2,...,  it  holds  that 


(s[ ,  t  >  0}  =4>  {—7 1  +  y/c  Bt,  t  >  0} 


where  {Bt,t  >  0}  is  a  standard  Brownian  motion. 

We  finally  come  to  the  bona  fide  long-range  dependent  case,  where  the  infinite 
sum  of  correlations  =  +00  is  induced  by  a  a  that  has  a  tail  behaving  like 

h(x)x~a,  1  <  a  <  2,  where  h(x)  is  a  slowly  varying  function.  As  a  result,  cr  belongs 
to  the  domain  of  attraction  of  a  non-Gaussian  stable  law  ([8],  [10],  [22])  with  index 
a,  a  fact  denoted  by  c t  EV (a). 
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Proposition  4.3  (Long-range  dependence)  Assume  that  a  €  T>(a),  1  <  a  <  2,  so 
that  P  [<j  >  x]  ~  h(x)x~a,  where  h(x)  is  slowly  varying.  By  taking  (r  to  be  the  1  fa- 
regularly  varying  function  satisfying 

r!i“  ^MCr)  =  (a-i)EM 


it  holds  that,  V  t  >  0, 

srt  =^r  —7 1  +  Xt 

where  Xt  is  a  stable  Sa  ^etr(2  —  a)  cos(7r  ^))“,l,o)  RV. 


We  sketch  the  proofs  of  Propositions  4.1,  4.2  and  4.3  in  Section  6,  where,  under 
long-range  dependence,  we  take  on  a  representative  Pareto-a  case.  A  short  discus¬ 
sion,  however,  of  Proposition  4.3  is  in  order  here.  It  is  well  known  (see  e.g.  [12], 
[2])  that  the  convergence  of  a  normalized  partial  sum  process,  such  as  >  0}, 
can  only  be  to  a  self-similar  process,  and  that  the  Hurst  parameter  H  may  be 
determined  through  the  regularly  varying  scaling  (r  by 


lim  =  xH ,  x  >  0 

r-voo 


Recall  that,  as  it  is  mentioned  in  [14]  and  [18],  with  a  Pareto-a  cr  the  M\GI\oo 
process  already  possesses  the  so-called  2nd  order  asymptotic  self-similarity  property, 
with  parameter  (3  —  a)/2.  That  is,  by  aggregating  the  original  process  {bn,n  = 
0, 1, . . .}  in  blocks  of  size  m  and  dividing  by  the  block  size,  one  can  obtain,  in  the 
limit  as  m  goes  to  infinity,  the  same  correlation  function  as  that  of  a  fractional 
Brownian  motion.  Because  of  this  fact  it  would  be  tempting  to  think  that,  except 
perhaps  for  a  slowly  varying  factor,  the  appropriate  scaling  squared  is  given  by  the 
rate  of  growth  of  the  partial  sums  variance,  r3~®,  and  that  convergence  occurs  to 
a  fractional  Brownian  motion.  Proposition  4.3  however  shows  that  this  is  not  true 
and  that  the  asymptotic  2nd  order  self-similarity  property  is  misleading:  The  proper 
scaling  does  not  exhibit  any  r^3-®)/2  dependence;  instead  it  contains  the  r1/®  factor 
associated  with  the  limit  theorem  for  the  RV  <r,  when  cr  €  T>(a),  where  1  <  a  <  2. 
The  marginal  distribution  of  the  limiting  process  is  not  Gaussian,  but  stable  with 
index  a,  1  <  a  <  2.  So,  the  limiting  process  {st°,t  >  0}  is  not  a  fractional  Brownian 
motion  but  an  a-stable  self-similar  process  with  infinite  variance,  Hurst  parameter 
H  =  1/a  and  dependent  increments,  i.e.,  it  is  not  the  a-stable  Levy  motion,  for 
which  the  increments  are  independent.  Further  work  on  determining  the  dependence 
structure  of  this  process  is  currently  in  progress. 
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5  Implications 


We  dealt  with  a  model  of  a  multiplexer  fed  by  M\GI\oo  processes,  that  account 
for  strong  dependencies  in  the  input  cell  stream.  Our  study  contributes  to  the 
analysis  of  such  non-standard  queueing  systems,  with  long-range  dependence,  in 
the  heavy  traffic  regime,  where  traditional  queueing  has  provided  useful  diffusion 
approximations.  Related  work  is  also  reported  in  [5],  [16]  and  [23].  Note  that  in 
the  work  of  Norros  ([16])  the  presence  of  fractional  Brownian  motion  is  postulated, 
while  in  [5]  Brichet  et  al.  show  how  fractional  Brownian  motion  can  arise  naturally 
in  the  heavy  traffic  analysis  of  on/off  sources  with  long-range  dependence.  Given 
that  the  M\GI\oo  processes  also  arise  as  the  superposition  of  infinitely  many  on/off 
sources  (see  [14]),  it  behoves  one  to  ask  whether  our  results,  that  belie  the  presence 
of  fractional  Brownian  motion,  contradict  those  of  [5].  This  apparent  contradiction 
is  resolved  though,  by  observing  that,  apart  from  the  different  way  in  which  heavy 
traffic  is  achieved  in  [5],  the  link  to  fractional  Brownian  motion  is  established  for 
a  particular  choice  of  the  on/off  periods,  yielding  a  subset  of  on/off  sources  that 
is  distinctly  different  from  that  considered  in  [14],  which  leads  to  the  M\GI\oo 
processes.  Thus,  our  results  stress  that,  in  contrast  with  short-range  dependence, 
within  long-range  dependence  itself  there  exist  modeling  possibilities  that  are  of 
fundamentally  different  nature.  Some  are  attracted  to  fractional  Brownian  motion, 
giving  rise  to  Weibullian  buffer  asymptotics,  yet  others  are  not;  we  have  shown  that 
the  M\GI\oo  model  is  attracted  to  a  a-stable  self-similar  process,  where  the  buffer 
asymptotics  are  expected  to  decay  like  x 1-“.  And  in  that  sense,  it  seems  that  the 
question  as  to  what  is  the  effect  of  long-range  dependence  does  not  have  a  clear-cut 
answer,  since  this  may  depend  heavily  on  the  particular  choice  of  model. 

6  Proofs 

We  provide  here  an  outline  of  the  calculations  that  determine  the  marginal  distri¬ 
bution  of  the  limiting  process  {s£°,i  >  0}.  First,  under  Assumption  3.1,  it  is  easy 
to  conclude  that 

lim  E [»r]  =  -7 1,  0  <t<T 

To  alleviate  notation  we  will  look  at  t  =  1;  any  other  t  can  also  be  handled  similarly. 
To  that  end,  let  us  consider  the  characteristic  function 

ipr(9)  =  E  [exp (*0Si)] ,  9  G  1R,  r  =  1, 2, ...  (6.1) 
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and  the  auxiliary  quantity 


Ar(0)  =  InE  exp  ^  (6^  —  ArE [cr])J  ,  0  g  R,  r  =  1,2,...  (6.2) 


By  writing 


4  =  r  E  (6»  -  A'E  H)  +  f (ArE  M  -  <0,  r  =  1, 2, . . .  (6.3) 

n=  1 


it  follows  that 


.(0)  =  exp  ^Ar(0)  +  30^-(ArE  [cr]  —  c)^  ,  r  =  1,2, .. . 


In  order  to  avoid  a  degenerate  limiting  RV,  we  are  interested  in  selecting  r  — 
1,2,...,  such  that  the  limit 

A(0)  =  lim  Ar(0),  0  e  1R  (6.5) 

r— >-oo 

exists,  is  non-zero  and  finite,  in  which  case,  under  Assumption  3.1,  so  does  the 
limiting  characteristic  function 

•0(0)  =  lim  i/v(0),  0  G  1R  (6.6) 

r— >oo 

So,  provided  that  the  limit  (6.5)  exists  and  Assumption  3.1  holds,  by  using  (6.4), 
ip (6)  can  be  expressed  as 

</,(0)  =  e-*>7+A(0),  Q  e  ]R  (6.7) 

After  some  algebra,  Proposition  2.2  enables  us  to  write 
A r(0)  =  —  ArE  [a]  +  ArE  [ct]  r  ^exp  —  1  j 

.  { 1  +  (j  _  exp  (_0  )  ±  (l  -  2)  r„  exp  (£»)  } 


=  —  ArE  [a]  r—  +  ArE  [cr]  i6 
Cr 


(«- 


f  1  -  exp  (-g)  *0  !  '  fie 

•U  + - i r — -  t  -E(r-n)r«exp(^n 

z  V  r  n=  1 


=  ArE  [a]  iO  — 
Cr 


(f)-1  A 


ArE  [a]  02 


(f)-1 


1  -  exP  ("?)  1  \ t,  (W  ^ 

- iT - -  72  2^(r-n)rnexpl  —  n\ 

Z  ’r  n=l  VC>r  ' 
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To  ease  the  computations,  we  introduce  the  quantities 

tfr(0)=r(l  +  £—  exP(0),  0€IR,  r  =  1, 2, . . .  (6.9) 

and 

Fr{0)  =  6>2-^  ^(r  -  «)rnexp  (j-n\  ,  9  €  1R,  r  =  1,2,...  (6.10) 

W  n=1  \sr  / 

In  all  cases,  we  will  seek  £r,  r  =  1,2,...,  ensuring  that  the  limits 


AT(0)  =  lim  Kr(9),  9  6  M 


(6.11) 


and 

F(0)  =rlhn  Fr(f),  6>GlR  (6.12) 

exist  and  are  finite.  Since,  by  Condition  3.1,  lim  (r  =  +oo,  it  follows  that 


lim 

r—yoo 


-1 


lim 

r— Kx> 


1  —  exp  (-ff) 
il 

Cr 


1 


(6.13) 


and,  by  virtue  of  Assumption  3.1  and  (6.9)  -  (6.12),  the  limit  (6.5)  of  (6.8)  can  be 
rewritten  more  compactly  as 


A(0)  =  -c  (K(0)  +  F{9)) ,  9  e  R  (6.14) 


As  we  have  already  mentioned,  the  choice  of  r  =  1, 2, . . .,  should  yield  K(9)  and 
F{9)  such  that  F(9)  +  K(9)  is  not  identically  zero. 

Short-range  dependence  Tw  <  +oo) 

With  the  scaling  ( r  =  y/r,  r  =  1,2,...,  (6.9)  and  (6.11)  yield 

if(e)  =  rKmr(|;-o(^))  =|.  (€K  (6.15) 

To  compute  F(9),  recall  first  that,  since  Tn  >  0,  n  =  1, 2, . . .,  the  implication 


OO 


E 


rn  <  -j-oo  =>  lim  nrn  =  0 

n—HX) 


holds  true  and,  by  Cesaro  convergence,  we  then  have 


lim  - 

r— >oo  t 


9  €  R 


(6.16) 


(6.17) 
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Applying  this  fact  to  the  limit  (6.12),  where  in  (6.10)  we  have  inserted  the  scaling 
Cr  =  Vr,  gives 

F(9)  =  92  ^lim^  53  exp  ^  €  1R  (6.18) 

In  order  to  show  that  the  exponential  factor  can  effectively  be  replaced  by  1  in  F(9) 
above,  we  form  the  difference  below  and  bound  it,  for  r£  =  1, 2, . . .,  by 


tr^exp^n)-1) 

<  t  r»  J  «p  ($«)  - 1 1  +  t- r»  |  «p  ($")  - 

JA  /  i9  \  " 

^  E  exp ( v9n) - 1  +  2Er« 

n=l  xv  '  n-re 


(6.19) 


Now,  since  for  fixed  rs  it  holds  that 

T‘ 

r-too  ■ 


lim  53  exp  ( —Ln)  —  1 


=  0,  9  e  1R 


by  invoking  the  fact  that 

00  oo 

53  Fn  <  +00  V  £  >  0 

n=l  n=rc 


for  re  large  enough,  we  see  that  each  of  the  two  terms  in  (6.19)  can  be  made  arbi¬ 
trarily  small,  so  as  to  yield 


lim 

r-¥  oo 


9  €  1R 


therefore  F(9)  of  (6.18)  is  simply 


OO 

F{9)  =  92  53  9  G  1R 

n=l 


(6.20) 


Substituting  (6.20)  and  (6.15)  in  (6.14)  and  inserting  this  in  (6.7)  we  obtain  a 
Gaussian  characteristic  function 


'ij)(9)  -  exp 


1  oo  \ 

— 02c(l  +  2  53  rn) )  >  -0eR 
Z  n=l  / 


(6.21) 
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Long-range  dependence  (X)^Li  Tn  =  +oo) 

We  illustrate  the  results  under  long-range  dependence  by  focusing  on  a  a  that 
follows  a  Pareto-like  distribution  with  parameter  a ,  where  1  <  a  <  2.  In  particular, 
we  take 

P 

for  which 

E  H  =  i  _2i-« 

and 

P  \o  >  nl  ,  r  , 

lim  — - - -  —  (a  —  1)E  [a] 

Consequently,  the  distribution  of  the  forward  recurrence  time  a  associated  with  a 
is  given  by 


[a  >  n]  = 


(n  +  I)1"*  -  (n  +  2)l~a 
1  -  21_“ 


n  =  0, 1, . . . 


r”spp>"]=Ep5;pi"  >  j)  =  («  +  l)1  a,  n  =  0,1,. 

t  J  =n. 


3—n 


We  begin  by  noting  that  the  case  a  =  2  needs  to  be  treated  separately  (Proposi¬ 
tion  4.2),  because  even  though  E  [a2]  is  infinite,  a  still  belongs  to  the  domain  of 
attraction  of  the  normal  distribution  (but  not  in  the  domain  of  the  so-called  normal 
attraction  of  the  normal  distribution,  for  which  E  [a2]  <  +oo)  .  The  input  pro¬ 
cess  {brn,n  =  1, 2,  •  •  •}  exhibits  long-range  dependence,  yet  is  barely  does  so,  since 
X)n=i  Pn  ~  lnr  =  h(r)  is  slowly  varying.  So,  undertaking  the  a  =  2  case  first,  we 
proceed  to  show  that  the  appropriate  heavy  traffic  scaling  is  —  y/rh(r)  =  VVlnr, 
r  =  1,2,  —  We  begin  by  substituting  this  scaling  in  (6.9),  so  that  (6.11)  reads 


K(0) 


lim  r 

r— ►oo 


id 

Vr  lnr 


=  lim  r 

r—¥  oo 


e2  Jo2  \\ 

2r  lnr  l  rlnr  I J 


lim  - — 

r— ►oo  lnr 


9  G  R 


(6.22) 


Next,  consider  (6.10),  which  in  this  case  becomes 
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«2  1  1  ^  (  i9n  \ 

-6Z- - O2— —  >  - -  exp  - 

lnr  rmr~ln  +  l  Vvrlnr/ 


for  r  =  2, 3, . . .,  0  €  IFt.  Since 
1 1 — l 


n  (  iOn  \ 

S^TTexpteJ 


<  r  -  1,  e  e  1R 


it  is  clear  that  only  the  first  sum  will  yield  a  non-zero  contribution  to  the  limit,  so 
we  have 


Fr(8) 


_2  1  "  1  (  idn  \ 

lnrUTin  Jo 


d/s/rlnr  T  \ 

Y  einx  dx 

71=1  ) 


lnr  Jo  elx  —  1 

j  /‘9/v/rinr  g irx  _  ^ 


dx 


(1  /•v/vrinr  pitx  _  i  \ 

1  +  i-r—  - - J-  dx),  9  GTR 

lnr  Jo  1  -  e~tx  J 


(6.23) 


One  can  verify  that  the  denominator  in  the  integral  above  can  be  rewritten  as 

1 

2 


1/(1  —  e~tx)  =  \  (1  —  icot  (x/2))  and  then  it  is  easy  to  see  that,  because  \e1TX  —  l|  < 
2,  we  get 


lim 

fO/Vr\nr  . 

/  ( enx  -  1 

)  cte 

r—¥  oo 

Jo  ' 

/ 

m 


<  lim  , - 

r-loo  yjr\ nr 


=  0,  9  e  1R 


so  that  Fr(9)  of  (6.23)  is 


Fr{9)  1  +  ^  (jrx  _  ^  cot  (a;/2)  dx^j  )  Q  e  R  (6 .24) 

In  Lemma  A.l  we  essentially  show  that  (6.24)  asymptotically  remains  unchanged  if 
we  replace  ^  cot  {x/2)  by  j,  thus  leading  us  to  the  expression 


F{9)  =  92  ( 1  +  lim  -i-  f 
y  r-+oolnryo 


j  rB/s/r lnr  gira:  _ 


dx 


=  02  1  +  lim 


1  eix  -  1 


r—too  lnr  Jo 


l 


X 


dx  ,  0  £  IR 


(6.25) 
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The  evaluation  of  the  limit  is  supplied  in  Lemma  A.2.  Inserting  the  result  in  (6.25), 
combining  with  (6.22)  in  (6.14)  and  finally  in  (6.7)  yields  once  more  a  Gaussian 
characteristic  function  : 

if)(6)  —  exp  ,  0  €  1R  (6.26) 

We  subsequently  study  the  case  where  1  <  a  <  2  (Proposition  4.3),  noting 
that  a  is  now  in  the  domain  of  attraction  of  a  non-normal  stable  law.  Of  primary 
importance  in  this  long-range  dependence  setup  is  the  parameter  H,  which  we  define 
as 

H  =  -  ,  1  <  a  <  2  (6.27) 

a 

We  will  show  that  the  appropriate  scaling  turns  out  to  be  (r  —  rH ,  r  =  1,2,..., 
giving  rise  to  a  non-Gaussian,  cr-stable,  self-similar  process  {s/°,  t  >  0},  whose 
Hurst  parameter  H  is  given  by  (6.27).  Observe  that  1  <  a  <  2  will  yield  H  in  the 
range  1/2  <  H  <  1,  while  in  our  previous  calculations,  where  a  =  2  ,  the  emerging 
self-similar  process  was  a  standard  Brownian  motion,  whose  corresponding  H  is 
exactly  1/2. 

Starting  from  K(6),  we  have 


K(0) 


Turning  to  Fr(Q),  r 
Fr(0)  = 


=  2, 3, . . .,  we  write 

-  n)(n  +  i)1-0  exP 

02  J2(r-n  +  l)nl~a  exp  (il(n  - 

n= 2  v ' 

°2^h  £(™1_a  -  n2"“)exp  f^(n  - 

n=l 

+e2^nr T  ^  E  nl~aexp  (JL  (n  _  ^ 


Since  1  <  a  <  2,  we  have  0  <  2 H  —  1  <  1,  therefore  the  third  term  above  vanishes 
as  r  — ¥  oo  and,  by  Cesaro  convergence,  so  does  the  second.  We  can  then  rearrange 
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the  only  contributing  summand  in  the  limit  (6.12)  as 


Clearly,  since  F(9)  and  F{—6)  are  complex  conjugates,  it  suffices  to  carry  out  the 
calculations  only  for  9  >  0.  Introducing  the  staircase  function  Tr  :  H  i->  IR+ 

Tr(x)  =  {  ^e<x^7^  0>0  n  =  1,2,... 

[  0  otherwise 

and  letting  kr  =  r1_ff,  r  —  1,2,...,  it  is  not  difficult  to  check  that  F(9 ),  0  >  0 

assumes  the  form 

F{9)  =  0afim  j\r{x)l-aexp{ikr9Tr{x))  dx  (6.29) 

^  Tr(x)2~a  exp(ikr9Tr(x))  dx,  9>  0 

Application  of  Lemma  A.3  to  the  each  of  the  two  limits  in  (6.29)  implies  that  we 
can  replace  Tr(x)  by  x,  and,  as  a  result,  F(0)  becomes 

F(0)  =  Km  r^(«-1)(2-«)  V")  exp  (ikr9x)  dx ,  9  >  0 

By  defining  the  quantities 

Ii-a(9)  =  0a  lim  rH(a-l)(,2-<x)  f  x1_“exp  (ikr9x)  dx,  9>  0  (6.30) 

r-too  Jq 

and 

h-a{9)  =  9a~l  lim  rH(a-1X2_“)  [  x2_aexp (ikr9x)  dx,  0>O  (6.31) 

r-nx>  j0 

we  can  express  F(9)  as 

F(9)  =  h_a{9)  -  I2-a(0),  9  >  0  (6.32) 

Appealing  to  the  results  of  Lemma  A.4,  where  we  carry  out  the  detailed  calculation 
for  li-o(0)  and  /2-«(0),  and  recalling  that  F(9)  and  F(—9 )  are  complex  conjugates, 
we  finally  obtain 

F(0)  =  |6>1°T(2  —  a)  exp  ^7rsgn(0)^  ^  ,  9  e  IR  (6.33) 
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Upon  inserting  (6.33)  and  (6.28)  in  (6.14),  then  this  in  (6.7)  and  rearranging,  we 
collect  the  limiting  characteristic  function 


ip(9)  —  exp  ^—iO'y  —  c|0|QT(2  —  a)  cos(n^  — — )  ^1  —  i  sgn(0)  tan(^)^  ^  ,  0eH 


(6.34) 


which  corresponds  to  a  stable  Sa  ^cF(2  —  a)  cos(7r^rp))  “  ,  1,  —7^  RV  (see  [22]). 


A  Appendix 


Lemma  A.l  If  lim  <5r  =  0  then 

J  r-+o o 


lim 

r— >  00 


eirx  -  1 
x 


dx  =  0 


(A.l) 


Proof.  The  Taylor  expansion  of  the  function  x  1-4  x  cot  x  gives  x  cot  x  =  1  —  x2/3  + 
o(x2),  which  enables  us  to  write 


1 

x 


( 


X  o(  X2)  (  1  o(x2)\ 

I2+V=*(i2  +  V) 


Now,  recall  that  Ve>03a;e>0  such  that  for  |a;|  <  x£  it  holds  that  o(x2)  / x 2  <  e 

and,  by  lim  5r  =  0,  3  re  such  that  when  r  >  re  we  have  |<ir|  <  xe.  As  a  result,  we 
r— >  00 

conclude  that 


1 

x 


X  X 

—  cot - 1 

2  2 


-  (ll  +  e)  |5r|’ 


r  >  re,  \x\  <  ST 


(A.2) 


Combining  (A.2)  with  \eirx  —  l|  <  2  and  with  our  assumption  lim  Sr  =  0,  quickly 

1  r— »-oo 

leads  us  to  (A.l).  ■ 


Lemma  A.2  It  holds  that 


lim  - — 

r— >00  inr 


0  ^  0 


(A.3) 
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Proof.  The  integrals  for  6  and  —6  are  complex  conjugates,  so  it  suffices  to  establish 
(A. 3)  for  6  >  0.  Let  us  define  R  =  and  write,  for  large  r, 


L 


R  eix  _  l 


dx 


R  eix  -  1 


x 


Jo  x  J 1 

/■l  pix  _  i  f-R  pix 

=  /  - - -dx-lnR+  /  — 

Jo  x  Jl  X 


dx 


(A.4) 


It  follows  that  only  the  In  R  term  above  will  offer  any  contribution  to  the  limit  in 
(A. 3).  To  see  this  note  that  the  first  term  is  easily  bounded  as 


fl  eix  -  1 
Jo  x 


dx 


<  1 


(A.5) 


while  the  last  one  can  be  handled  by  a  contour  integration  : 


/-l  pZ  pZ  rm  pz  rRein  pz 

-dz+  -dz+  -dz+  — 

-R  z  Jeilr  Z  Ji  Z  jRe »»/*  Z 

/—I  g*  i*7T /2  rR  g ix 

—  dx  +  i  /  exp  (cos  <p  +  i  sin  <j>)  d</>+  /  — 

-R  X  J-k  Jl  X 

/*7T 

+i  /  exp(f?(coS(£  +  ?sin</>))  dip 
J-nll 


0  = 


dz 

R  cfryj 


dx 


This  implies  that 

rR  pin 


rR  pi#  rR  p—%  fit  .  fit 

/  —dx  <  /  —dx+  ecos  <t>d(f)+  i 
Jl  X  Jl  X  Jtx/ 2  Jir/2 


0R  cos  </> 


d(f) 


(A.6) 


where  all  the  terms  in  the  right  hand  side  are  bounded  uniformly  in  R.  Then,  because 
of  (A.5)  and  (A.6),  (A.4)  will  yield  the  desired  conclusion  (A. 3),  after  multiplying 
by  1/lnr,  taking  the  limit  r  ~¥  +oo  and  using  the  fact  that 


i"K/iS)  i 

im  - - — - - —  — - 


lim 

r-+oo 


lnr 


2’ 


9  >  0 


Lemma  A. 3  For  1  <  a  <  2,  p  >  — 1,  H  =  ^  and  6  >  0  it  holds  that 


lim  rtf(«-i)(2-«) 
r— t-oo 


/  (Tr(x)p  exp  ( ikr6Tr(x )) 
Jo 


—  xp  exp  ( ikr6x ))  dx  =  0 


(A.7) 
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Proof.  We  first  bound  the  difference  by 

[0  rd 

/  Tr(x)p  exp  (ikr6Tr(x))  dx  —  I  xp  exp  (ikT0x)  dx 

Jo  Jo 

I  [9 

<  \  xp  (exp  (ikr9Tr(x))  —  exp  (ikrBx))  dx 

\Jo 

re 

+  /  (Tr(x)p  —  xp)  exp  (ikr6Tr(x))  dx  (A.8) 

Jo 

and  then  bound  each  of  the  two  terms  separately.  Recalling  that  kr  =  rl~H  and 
Tr(x )  —  x  <  r  =  1,2, . . the  first  term  is  bounded  as  follows  : 

[9 

/  xp  (exp  (ikr6Tr(x))  —  exp  (ikr0x))  dx 

Jo 

fe 

<  xp  |exp  ( ikr6(Tr(x )  —  x))  —  1|  dx 
Jo 

f0  00  un 

<  /  xp  'V'  —  | Tr(x)  —  x\n  dx 

-Jo  t,n! 

r0  00  t.n  on 

<  xp'V'  dr 

Jo  n!  rn 

=  —  l)  [  xp  dx 

“  £  (r:~  0  (A.9) 

We  multiply  by  r#(a-1)(2_“)  and  calculate  the  limit  : 

^lim)r^a_1^2_a^e0r  H  —  l)  =  lim  rH^a~l^‘2~c^  {0r~H  +  o  (r~^)) 

rH(a- l)(2-a)  f  0  ( r~H )  \ 

=  Hm  - jj -  6  +  -A-„  > 

r—¥oo  qrn  1  rp—fi  I 


=  lim  r  H(2  ^2)  \0  + 


o  (r 


(A.10) 


For  the  second  term  in  (A.8),  with  p  >  —  1,  we  have 


f 9  (Tr{x)p  -  xp)  exp  (ikr6Tr(x))  dx  <  [ *  {Tr(x)p  -  xp)  dx 

Jo  Jo 
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< 


< 


m-r- 
l 

)P+l  (r  +  l)p+1 
^+T  p  + 1 

9P+ 1 


=1 

9p+i 

rp+i 


xp  dx  — 


0P+\ 
p  +  1 

QP+l 


+  1 


-r 


p + 1 

QP+l 

P  +  1 
1 


(A. 11) 


Again,  multiplying  by  rH^  x)(2  “)  and  writing  (1  +  y)p+1  -  1  =  +  o  we 

calculate  the  limit  as 


lim  rH(«-l)(2-a) 

) — >oo 


P+1 


lim  r#(®-l)(2-a) 
r— ►  oo 


„H(a— 1)(2— a) 

lim  - - 

r-*-oo  f 

0 


(A.12) 


where  the  last  equality  is  derived  easily  by  observing  that,  whenever  1  <  a  <  2, 
rH{a-i){2-a)  _  After  multiplying  (A. 8)  by  r^(“-!)(2-a)  an(j  taking  the  limit 

as  r  ~+  +oo  we  invoke  (A.12),  (A. 11)  and  (A. 10),  (A. 9)  to  obtain  (A. 7).  ■ 


Lemma  A. 4  For  1  <  a  <  2  it  holds  that 

Ti-a(0)  =  9aT(2  —  a)  exp  >  0 

and 

h-a(0)  =0,  0  >  0 


(A.13) 

(A.14) 


Proof.  Both  statements  are  clearly  true  when  0  =  0.  To  prove  (A.13)  for  6  >  0 
we  start  from  definition  (6.30)  and  perform  a  change  of  variable  u  =  krx  to  get 


T0)  =  9a  lim 

r-kx> 


rOkr 

J  0 


ul~aeiud  du, 


6  >  0 


(A. 15) 
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Writing  R  =  9kr,  r  =  1, 2, . . for  the  upper  limit  and  observing  that  lim  R  =  +oo, 
0  >  0,  this  integral  can  be  evaluated  by  the  standard  first  quadrant  contour  inte¬ 
gration  below: 


/•R  .  i-Re1*12  riro  /*roe*° 

=  /  zl~aelz  dz+  z1^aetz  dz  +  zl~aetz  dz+  z 1 

Jr0  J  Rei0  JiR  Jr0e *"/* 

=  J  xl~aetx  dx  +  iR2~a  J  exp  (i<p{2  —  a)  4-  iRe1^  dp 

rro  /»Q  . 

+i 2~a  J  yl~ae~y  dy  +  irl~a  J  exp  (i<p{2  —  a)  +  J  dp 


V*  dz 


Letting  ro  go  to  0,  we  have 
rR 


J  x1  aetx  dx  =  —iR2  a  J  exp  (ip( 2  —  a)  +  iRe1^  d(f> 

+i2~a  [  yl~ae~y  dy 
Jo 

Observe  that  the  last  integral  above  yields  a  Gamma  function  : 

i2~a  lim  [  yl~ae~y  dy  =  —e~l~Y(2  —  a),  6  >  0 

r-»oo  J0 


(A.16) 


(A.17) 


Next,  calculate  the  limit  involving  the  other  remaining  integral  in  (A.16),  using  the 
bound 


(ip(2  —  a)  +  iRe1^  dip 


< 


exp  (— Rsinp)  dcp 


(A.18) 


and  the  fact  that  sin^>  >  0  <  <p  <  7r/2,  which  is  implied  by  the  concavity  of 

py-^smp-^-,  for  <p  €  [0,  f]  : 


rir/2 

lim  R2~a  /  exp(— Rsinp)  d<p  < 

r-+oo  Jo 


lim  R2~a 

r—too 


-  lim  R2~a 

2  r-±oo 


—  lim  Rl-a  (l  -  e~R) 
20  r-t  oo  V  J 

o,  e  >  o 


dcp 


(A.19) 


Taking  the  limit  as  r  — ►  oo  of  (A.16),  combining  (A.17),  (A.18),  (A.19)  and  inserting 
in  (A.  15)  readily  yields  (A.  13)  for  6  >  0. 
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Having  shown  (A.  13),  it  is  easy  to  prove  the  claim  for  /2_a(0)  by  integrating  by 
parts  : 

rd 


f* 


,2 —OlgikrX 


dx  =  —  02"VM 

1  rCf*  y 


(2  -  a)  xl~aeikrX  d^j  (A.20) 


After  multiplying  the  equation  above  by  rH(a-i)(2-a) ,  (A.  14)  follows  from 

7  rH(a- l)(2-a)  1 


\h-a(0)\  <  ea~l  lim 

l  r-too  kr 


04~a+  lim  —  (2  -  a)0-a\h-a(9)\ 
*■->  oo  kr 


=  0  lim  r  i2TJ-  +  1(2  -  a)\h^a(0)\  lim 

r-¥oo  n  v  ’ 1  r-i-oo  kr 

=  0,  6  >  0 
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